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Abstract 

The particle density distribution emerging from the solution of the Vlasov equation 
for relativistic, non-interacting particles with spherically symmetric initial conditions is 
shown to exhibit a shell-like structure for late fixed times in the center-of-mass (CM) 
frame of the system. A similar phenomenon was recently observed employing the test- 
particle method to solve the Vlasov equation for quarks with Nambu-Jona-Lasinio-type 
interactions, and was attributed to the attractive forces among the particles. Contrary 
to this claim, it is demonstrated here that this effect is of purely relativistic origin and 
is sensitive only to the mass of the particles and their initial phase-space distribution. 
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One of the primary goals in present relativistic heavy-ion collision experiments at the AGS 
at Brookhaven, the SPS at CERN, and in future experiments at Brookhaven's Relativis- 
tic Heavy Ion Collider (RHIC) and CERN's Large Hadron Collider (LHC) is the temporary 
formation and subsequent observation of the so-called quark-gluon plasma (QGP), the decon- 
fined, chirally restored phase of strongly interacting matter [ffl, In the common picture of 
chiral symmetry restoration 0], at sufficiently high temperatures or densities the massive and 
confined quasi-particle excitations, the constituent quarks, become bare, undressed quarks 
with current quark masses being much smaller than the constituent quark masses. 

Recently, Abada and Aichelin || studied the dynamical evolution of a collision-free quark- 
antiquark system by solving the relativistic Vlasov equation for the one-particle distribution 
function /(x, p, t), 

% + V^-V p / = , (1) 

where the effective momentum-dependent force 

-V X E = -cV.[p 2 + M c 2 (x,t)c 2 ] 1/2 = ~ Mc(X F t)c2 c 2 V g M c (x,t) (2) 

is generated by a self-consistent gap equation for the constituent quark mass M c derived within 
the chiral quark model of Nambu and Jona-Lasinio j|, [|. The initial condition was a sphere 
of radius R, homogeneously filled with quarks and antiquarks having a thermal momentum 
distribution corresponding to a temperature T well above the transition temperature T c ~ 140 
MeV. Such "hot spots", i.e., small-scale fluctuations with sufficiently high (and thermalized) 
energy density, are expected to occur in present heavy-ion experiments at the SPS |J and 
also in future collider experiments at RHIC @. 

The striking observation of Abada and Aichelin was that, in the evolution of the system, 
most particles remain concentrated in an expanding shell, and that consequently the phase 
transition takes place at the inner and outer surface of that shell. Thus, the situation is fun- 
damentally different from the naive expectation that the initial fireball expands by retaining 
its shape and consequently hadronizes from its surface only. For the interpretation of the re- 
sults, Abada and Aichelin argued as follows: in the non-relativistic case the thermal velocity 
distribution of particles peaks around (v) = ^2T/m, and in an interaction-free expansion, 
all the many particles with that velocity will stay together, leading to a corresponding peak 
in the density, i.e., the observed shell-like structure. For interacting quarks, the attractive 
forces in the Nambu- Jona-Lasinio model will enhance this effect by trying to keep the (bare) 
quarks in the region of high density. This conclusion was further motivated by a similar 
observation found in non-relativistic simulations of expanding nuclear matter, where a shell 
structure forms after the nucleons have frozen out and enter the liquid-gas transition regime 

In contrast, in this note we show that the observation by Abada and Aichelin ||, i.e., 
the shell-like structure in the density distribution, originates simply from ultrarelativistic 
kinematics. Moreover, contrary to their line of arguments, we show that in the non-relativistic 
limit, the density distribution does not exhibit such a structure. We will not justify the 
validity of the Vlasov equation for modelling the expansion of a quark-antiquark plasma. 
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Obviously, collisions among the constituents are important to achieve the assumed, thermally 
equilibrated initial condition in the first place. The relevant quantity for thermalization is 
given by the transport rate r* rans of the quarks in a QGP 0. Estimates with a simple 
infrared cutoff suggest a rather high rate of r* rans « 0.3 — 0.5 T (T denotes the temperature 
of the plasma), however, a more detailed treatment of the infrared sector in those rates by 
means of resummation methods raises the question whether one can trust those perturbative 
calculations at all ||. 

The collision rate may be high enough to maintain (local) thermal equilibrium even during 
the expansion. In that case the system's evolution is governed by ideal hydrodynamics. 
We shall compare the relativistic hydrodynamical solution for an expanding sphere with the 
results obtained from solving the Vlasov equation. 

On the other hand, it was also recently speculated that, in particular shortly before the 
phase transition happens, the collisions in the plasma become so rare that their effect on the 
overall dynamics can be neglected. Then, the system freezes out already before the transition 
occurs and it is appropriate to model the non-thermal dynamics with the Vlasov equation 

There is a very simple, intuitive explanation why a shell occurs in the relativistic expansion 
of the fireball: expanding matter piles up in a shell, because the velocity of light c limits 
the maximum particle velocity. Suppose the particles are ultrarelativistic, i.e., massless. In 
that case, the velocity of all particles equals the velocity of light, v = c. Obviously, without 
collisions, at times t ^> R/c they can only be located in a spherical shell of thickness 2R, 
moving outward with light velocity. (The particle density in that shell is not necessarily 
constant, but depends on the initial distribution of momenta and coordinates in the sphere, 
see below). The interior of the shell will be completely void of particles. For particles with 
finite mass, the thickness of the shell is roughly 

d « 2R + Av t , 

where Av = \(v 2 ) — (v) 2 ] 1 ^ 2 characterizes the average spread of the velocity distribution. 
That spread grows with increasing mass of the particles, such that in the non-relativistic 
limit the interior of the shell is filled with particles, thus restoring the naive expectation for 
the expansion of a sphere (and disproving the arguments presented in H). 

As will be argued below, the modification of the momenta and velocities of the particles 
due to the force term (|2|) and the time dependence of the mass are small, so that the generic 
features of phase-space dynamics will be represented by the free-streaming limit of eq. ([[]), 

f+v.V,./ = 0, (3, 

where v = pc 2 /E and the mass M c of the quarks is assumed to be independent of time. In 
the following we explicitly construct the solution for the particle density distributionF] 

pCM) = /(^/( X 'P>*)- ( 4 ) 

1 We choose units Ti = ks = 1. 
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Our initial condition will ultimately be the same as in ||, i.e., the particles are thermally 
distributed in momentum space and homogeneously distributed in coordinate space within a 
sphere of radius R. For the moment, however, arbitrary (though homogeneously filled) initial 
coordinate-space volumina V and arbitrary initial momentum distributions n(p) are allowed. 

We introduce the characteristic function of the initial volume, Oy(x) = 1 for x G V, and 
elsewhere. Then, for a coordinate-space distribution which is homogeneous in V, the initial 
phase-space distribution separates as 

/(x, P ,o) = e v (x)«( P ) 

= y>ae y (a)/ a (x,p,0) (5) 

where we defined the phase-space distribution for a 5-like source, 

/ a (x,p,0) = 5 3 (x-a)n( P ) . (6) 
Since ([|) is a homogeneous linear differential equation, the solution /(x, p, t) has the form 

/(x,p,t) = |d 3 ae v (a)/ a (x,p,t) . (7) 

The solutions f a to eq. ([3|) (for time-independent particle mass) are readily obtained: 

/ a (x, p, t) = / a (x - vt, p, 0) = 5 3 (x - vt - a) ra(p) . (8) 

The particle density distribution (£|) p a (x, t) for a 5 like source @ located at x = a follows 
after a careful evaluation of the 5-function as 



p a (x,t) = / ^—^ n(p) 5 3 I x - — i - a 

1 , . s , s /m x 3 



n(m cosh 77 u) cosh r; — I 0(1 — tanhr/) , t>0, (9) 



where 



(2tt) 3 v ' ' '\t 

u ee , tanhr; ee \u\/c . (10) 

t 

Obviously, m cosh r\ u is the momentum and mc 2 cosh 77 the energy of the particles. The step 
function assures causal expansion. The non-relativistic limit corresponds to |u| c, cosh 77 w 
1 and becomes 

Pa r -(x,t) = ^ F 77(7nu)(^) 3 . (11) 

If the initial momentum distribution n(p) is isotropic, i.e., n(|p|), and a monotonously decreas- 
ing function of |p|, as it is the case for a thermal (Bose-Einstein, Fermi-Dirac, or Maxwell- 
Boltzmann) distribution, the density distribution p a (x, t) is also isotropic, p a (r,t), and, in the 
non-relativistic limit, a decreasing function of the (radial) distance r ee |x— a| from the source. 
This disproves the argument presented in 0] about the formation of a shell-like structure in 
the non-relativistic limit. 



3 



This behaviour changes completely in the ultrarelativistic regime due to the additional 
factor cosh 5 rj. This term is a monotonously increasing function of r and, in combination 
with the momentum distribution function, thus leads to a pronounced peak in the density 
distribution. In Fig. 1 the quantity (ct) 3 p a (r,t) is plotted as a function of r/ct for various 
masses m. (Note that this representation is time-invariant.) A Fermi distribution 



n(|p|) = n F (\p\) = ^ E/T - - , E = c^p 2 + (mc) 2 , 

was chosen with a temperature T = 160 MeV, slightly larger than the expected critical 
temperature T c . The shape of the distribution is very sensitive to the ratio mc 2 /T. One 
clearly sees how the peak structure emerges already when mc 2 /T is of order one, as most of 
the particles are already moving with a velocity close to c. (If mc 2 /T <C 1 the width of the 
velocity distribution is proportional to ~ (mc 2 /T) 2 .) 

For a uniformly filled sphere of radius R, 

„ / s ( I for lal < R 
6 ^ (a) = { for|a|>i? • 

The particle density distribution p(x, t) is then obtained as 

p(x,t) = J d 3 a e v (a)p a (x,t) . (12) 
In principle, one should multiply (12) by the degeneracy d q = 12 of light quarks (and similarly 



for the antiquarks), but this factor is, for the sake of simplicity, omitted in the following. Eq. 
(|12"D can be reduced to a one-dimensional integral and takes the form 

(TTIC \ 3 7T ~ rmin(r+l;t) . . 

— ~) — Q(t — f + 1) / dzz ( 1 — (z — f) 2 ) n(mcsmhr] z ) cosh 5 r\ z 

2-KtJ r Jf-i v > 



for r > 1 ; 



(13) 



p(f,t) 



( mc\ 3 


(•min(l 


{—-) 2tt 


2 / 


\27TtJ 


Jo 



dzz n(mcsmhr] z ) cosh r] z 

| /•min(l+f;i) 



1 _ rmm(L+r:t) 

— Q(t + f— 1) / dz z ( 1 — (z — r) ) n(mcsinh?7 2 ) cosh 

2f Ji-f v ' 



for r < 1 , 

where the dimensionless variables r = r/R and t = ct/R have been introduced and tanh^ = 
z/t. Expression (|1^) has to be integrated numerically. 

In Fig. 2 the resulting density profile p(f,t) is shown for different times t between 0.1 and 
3. Again, the initial momentum distribution is a Fermi distribution np with a temperature 
T = 160 MeV. The mass of the quarks is taken as m(= M c ) = 50MeV/c 2 . For times t < 1 
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the density in the central region r < 1 — t remains constant due to causality. In other words, 
it remains constant because as many particles leave an infinitesimal volume element as do 
enter. However, since most particles move with nearly light velocity (the thermal velocity is 
(v) m 0.9858 c for the present situation, while the spread of the velocity distribution is only 
Av = \(v 2 ) - (v) 2 \ 1 / 2 « 0.0031c), for times * > 1 no particles can enter the central region 
around f = any more, so that it becomes completely depleted. Thus, a shell starts to move 
outwards at time t — 1, as one clearly recognizes in Fig. 2. As expected, that shell has a 
width of 2 R and travels with the velocity of light. 

In Fig. 3 density profiles p(f, 1.5) are shown for different masses. For all masses smaller 
than 300 MeV/c 2 relativistic effects become important, and a shell structure emerges in the 
temporal evolution. 

This behaviour is quite similar to the relativistic, ideal hydrodynamical expansion of a 



spherical fireball [ll|. Also in this case, the particle density distribution exhibits a shell 
structure at fixed times in the CM frame of the fireball. The reason is that moving matter at 
finite r experiences relativistic time dilation ]I2] with respect to matter in the center (which is 



at rest due to symmetry) and thus dilutes less rapidly. To confirm this, and to investigate the 
effect of a varying particle mass on the expansion, we solved the equations of ideal relativistic 
hydrodynamics, 

= (14) 

(T^ v = 4- p')u fl u u — pg^ v is the energy-momentum tensor of an ideal fluid, e, p are energy 
density and pressure in the local rest frame of a fluid element, moving with 4-velocity u M in 
the CM frame of the fireball, g^ u = diag(+, — , — , — ) is the metric tensor), for a spherically 
symmetric initial fireball of radius R and temperature T = 160 MeV and for ideal Fermi gas 
equations of state 

1 f°° i 

e(T) = ~~ ~ / dpp 2 \ lp 2 + m 2 c 2 n F (p) , (15) 
2tt z Jo v 

1 r°° p A 

P( T ) = T~2 / d P I 2 i 2 =7 ' ( 16 ) 

2ti z Jo VP + m c 

with masses m varying between 10 and 940 MeV/c 2 . The corresponding particle densities are 

1 f°° 

P(T) = ^2 J Q dpp2 Uf ^ ■ ( 17 ) 
The hydrodynamical equations ( |T3D are solved with the (one-dimensional) relativistic Harten- 



Lax-van Leer-Einfeldt algorithm fL2|, supplemented with a Sod predictor-corrector step 13 



to account for spherical geometry [14]. 

In Fig. 4 we show the corresponding density profiles at time t = 2.97. The situation 
is qualitatively very similar to Fig. 3, which is remarkable considering the fact that ideal 
hydrodynamics corresponds to the limit of an infinite collision rate (and thus immediate local 
thermodynamical equilibrium). The main quantitative differences are (a) the time t has to 
be about twice as large as before to obtain density distributions similar to Fig. 3, and (b) the 
quark densities are smaller. The reason for (a) is that for a time t — 1.5 the hydrodynamical 
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rarefaction has not yet reached the center of the fireball (the rarefaction velocity is the velocity 
of sound c s < c/v3 ~ 0.5774c, in contrast to the free-streaming case where the rarefaction 
wave travels with the velocity of light c), and thus the shell structure could not possibly have 
formed. The reason for (b) is that the quark density is not conserved in the solution of (0), 
in contrast to the conservation of the quark number in the free-streaming scenario: due to 
the work performed in the hydrodynamical expansion the system cools, and thus temperature 
and, consequently, particle density decrease. Thus, the densities are smaller in Fig. 4 than in 
Fig. 3. In order to conserve the quark number in the hydrodynamical evolution, one would 
have to supplement the equations of motion (|T4| ) with the continuity equation for the quark 
number current, but this would also require introducing a finite quark chemical potential in 
the equation of state. This is certainly beyond the scope of our present work, where we are 
interested in qualitative similarities between relativistic free-streaming and hydrodynamical 
solutions. 

Finally, we address the potential influence of the self-consistent mass term and the corre- 
sponding force (0) on the expansion dynamics. The mass M c (r, t) influences the velocities 
of the particles in two ways. First, due to its spatial dependence a force is generated which 
changes the momentum of a particle in time according to 

M c (r,t)c 2 9 „ „ , , 
p = pe p + pe p = ^c'V.M^r.i) . (18) 

Second, the velocity of a particle changes because of the time dependence of the mass; in 
particular it will continuously decrease if M c (t) increases with time as in the present case. 

According to || the initial mass of the particles is M c ps 50 MeV/c 2 . For a temperature 
of T ps 240 MeV (used in ||) the thermal velocity is ~ 0.9931c. To estimate the typical 
change in momentum, the modulus of the gradient in the mass M c is (conservatively) read 
off the results of || as |Va;M c (r, t) c 2 \ < 500 MeV/fm, acting over a period in time At ps 0.5 
fm/c. Taking (E) and (pc) ~ 3T and an average mass (M c c 2 ) ~ 100 MeV the change in the 
(average) momentum is estimated as 

(Ae p ) = (ep)At < \V x M c (r, t)c 2 \ At w 0.05 . 

Hence, the typical momentum of a particle is to good approximation conserved both in mag- 
nitude and direction; the change is at most a few percent. 

For approximate momentum conservation, pi(t = tj)c ps 3T ps pf(t = tf)c, the change in 
the typical velocity of the particles due to the increase of the mass from M c (ti) ~ 50 MeV/c 2 
to M c (t f ) pa 300 MeV/c 2 is readily obtained as 

^;Af c =50MeV/c 2 - f/ ; M c =300MeV/c 2 ~ 0.09 C . (20) 

All particles with large (thermal) velocities will be decelerated by typically ten percent. As 
a consequence, the formation of the shell will be somewhat delayed. Furthermore, this de- 
celeration of the particles will stabilize the region of expanding matter, where the mass of 
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the particles is still rather low, i.e., the interior of the outgoing shell. In conclusion, taking a 
self-consistent mass into account will only have minor effects on the dynamics. The generic 
structure of phase-space dynamics is dominated by relativistic kinematics as given by the 
solution of (|]). 

To conclude, in this note we have investigated the free-streaming expansion of a fireball 
consisting of quarks and/or antiquarks. We explicitly constructed the solution for spher- 
ically symmetric initial conditions and confirmed the existence of a shell-like structure in 
the expansion at late times [[|. Contrary to the arguments of [[|, we have explained that 
structure as arising solely from (and only for) relativistic kinematics. Interestingly enough, 
similar structures occur in the ideal hydrodynamical limit (where interaction rates are infi- 
nite), which again confirms that their origin is kinematical rather than due to specific details 
of the interaction between the constituent particles. 

The overall decrease in the density of the shell is ~ 1/t 3 , cf. eq. (0). If the critical density 
for the phase transition is reached when the shell was already formed (for times t > R/c, 
cf. Fig. 2), the transition necessarily happens at the outer as well as the inner surface of the 
shell. This is what has been observed in [[3]. In principle, however, it is also conceivable that 
the phase transition happens before the shell structure has formed, for instance when the 
initial temperature is only slightly above the critical temperature T c . Then the fireball will 
hadronize only in the surface region where the density decreases fast at times t < R/c. 

The analytic solutions presented here may also serve as benchmark tests for numerical al- 
gorithms like the test-particle method. Another test case with a solution that can be readily 
obtained by similar means (but will not be presented here for the sake of brevity) is the 
free-streaming expansion of a slab (finite thickness in one dimension, but infinite extension in 
the other two). In that case, no shell structure emerges (because sufficiently many particles 
enter the central region from other parts of the slab). This is in contrast to the analogous 



hydrodynamical expansion problem fT2|1 
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Figure captions: 



Figure 1: 

The density profile (ct) 3 p a (r, t) of a point source for various masses m as a function of r/ct. 
The initial momentum distribution is a Fermi distribution with temperature T = 160 MeV. 

Figure 2: 

The density profile p(f,t) at various stages t of the free-streaming expansion of a fireball. 
The initial momentum distribution is a Fermi distribution with T = 160 MeV. The mass of 
the particles is taken as m = 50 MeV/c 2 . At times t > 1 a shell structure emerges in the 
radial expansion. 

Figure 3: 

The density profile p(r, t) at a fixed time t = 1.5 for the free-streaming expansion of a fireball. 
The mass of the particles is varied between 10 and 940 MeV/c 2 to demonstrate the relativistic 
effects in the expansion. 

Figure 4: 

The density profile p(f, t) at a fixed time t = 2.97 for the ideal relativistic hydrodynamical 
expansion of a fireball. The equation of state is that of an ideal gas of relativistic particles 
with mass varying between 10 and 940 MeV/c 2 . 
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